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ABSTRACT 

The precision of intensity measurements of the extragalactic X-ray Background (XRB) 
on an angular scale of about a degree is dominated by spatial fluctuations caused by 
source confusion noise. X-ray source counts at the flux level responsible for these fluc- 
tuations, 10~^^ergcm~^ s""'^, will soon be accurately measured by new missions and 
it will then be possible to detect the weaker fluctuations caused by the clustering of 
the fainter, more distant sources which produce the bulk of the XRB. We show here 
that measurements of these excess fluctuations at the level of (^) ~ 2 x 10~^ are 
within reach, improving by an order of magnitude on present upper limits. Since it is 
likely that most (if not all) of the XRB will be resolved into sources by AXAF, sub- 
sequent optical identification of these sources will reveal the X-ray volume emissivity 
in the Universe as a function of redshift. With these ingredients, all-sky observations 
of the XRB can be used to measure the power spectrum of the density fluctuations in 
the Universe at comoving wavevectors kc ~ O.Of — O.lMpc^^ at redshifts where most 
of the XRB is likely to originate {z ^ 1 — 2) with a sensitivity similar to, or better 
than, the predictions from large-scale structure theories. A relatively simple X-ray 
experiment, carried out by a large-area proportional counter with a 0.5 — 2deg^ col- 
limated field-of-view scanning the whole sky a few times, would be able to determine 
the power spectrum of the density fluctuations near its expected peak in wavevector 
with an accuracy better than 10 per cent. 

Key words: Methods: statistical - diffuse radiation - large-scale structure of Universe 
- X-rays: general 



1 INTRODUCTION 

1.1 The X-ray Background and the Large-Scale 
Structure of the Universe 

Thirty-five years of study after its discovery by Giacconi et 
al (1962), the X-ray Background (XRB) has proven to be a 
valuable tool in the study of the Large-Scale Structure (LSS) 
of the Universe. A significant fraction of the XRB (particu- 
larly at soft energies) is now resolved into sources, the vast 
majority of which are extragalactic. The deepest surveys car- 
ried out with ROSAT (see, e.g., Hasinger 1996 for a recent 
review) resolved about 60 per cent of the 0.5 — 2keV XRB 
into mostly Active Galactic Nuclei (AGN) and other X-ray 
luminous galaxies (particularly Narrow-Line X-ray Emitting 
Galaxies-NLXGs) . The X-ray volume emissivity of these ob- 
jects rises rapidly from redshift z = 0to2 = l — 2 above 
which it appears to decline. It is at redshifts of 1-2 where the 
bulk of the resolved soft XRB originates. Preliminary work 



on ASCA observations (Georgantopoulos et al 1997) hints 
that something similar is happening at harder energies, al- 
though the fraction of the XRB resolved is lower and the 
sources much more difficult to identify due to the limited 
spatial resolution of ASCA. 

Isotropy arguments at soft X-ray energies (see, e.g., 
Carrera, Fabian & Barcons 1997) show that a good deal 
of the unresolved fraction of the soft XRB has to come from 
redshifts z > 1. It is then concluded that the bulk of the 
XRB originates precisely at the epoch where the largest 
present-day structures collapse and where different cosmo- 
logical models give the most different predictions. 

Any cosmological model has to confront two basic 
boundary conditions: (1) the Universe was very smooth at 
2 — 1500 where the Cosmic Microwave Background (CMB) 
was produced; and (2) the Universe is very lumpy and 
rich in structure today {z — 0). Density fluctuations in 
the Universe grow linearly from the highly homogeneous 
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CMB recombination epoch until they become non-hnear at 
low redshifts and collapse to form the LSS we see today. 
The redshift at which a density fluctuation becomes non- 
linear depends on the spatial scale of the fluctuation and 
on the cosmological model, but for the most popular mod- 
els, this happens at a typical redshift z ~ 1 — 5 for scales 
of 10 - 100/i^^Mpc {Hq = lOO/ikms^i Mpc-\ go = 0.5 
and A = are used throughout unless explicitly stated). 
Studies of the XRB as proposed here will add a further con- 
straint to cosmological models, since they will measure the 
power spectrum of the fluctuations at comoving wavevec- 
tors kc ~ 0.01 — 0.1/iMpc~^ with sufficient sensitivity to 
constrain cosmological scenarios. In particular the parame- 
ters governing the evolution of the power spectrum (e.g., qo) 
could be constrained by these observations. 

Mapping the structure of the Universe at intermediate 
redshifts with X-rays has the added bonus that the objects 
presumably occur in the highest density regions. Indeed, 
cosmic X-ray sources require strong gravity to be switched- 
on (either very extended potential wells, as in galaxy clus- 
ters, or very deep ones as in AGN), and therefore the X-ray 
sources existing at intermediate redshifts could actually cor- 
respond to the objects that first formed in the Universe at 
even earlier epochs. A possible indication of this is in the 
'bias factor' for X-ray sources for which there is some evi- 
dence that it has a high value, at least for the nearby lu- 
minous X-ray sources (see, e.g., Miyaji 1995 and Boughn, 
Crittenden & Turok 1997 who find values larger than 5). 
Although direct studies of the clustering of X-ray selected 
AGN (Boyle & Mo 1993, Carrera et al 1997) show that these 
values might be too large for lower-luminosity objects, bias 
factors of 2 or 3 might be likely for X-ray sources. This is- 
sue will be solved with ABRIXAS when the dipole of the 
distribution of sources is compared to the XRB dipole. 

Parallel studies at other wavelengths are currently be- 
ing used to extract information about LSS at high redshift. 
Deep galaxy surveys, and especially the Canada-France Red- 
shift Survey, are showing how the galaxy-galaxy correlation 
function evolves up to redshifts close to 1 (Le Fevre et al 

1996) . However, scales of the order ~ 100/i~"'^Mpc are be- 
yond the scope of such work. High redshift clustering is also 
being investigated by means of different classes of QSO ab- 
sorption systems (Fernandez-Soto et al 1996, Cristiani et al 

1997) . None of these studies provides clear evidence on how 
clustering in the Universe evolves. Again, scales in excess of 
10/i-^Mpc are not easily accessible by these studies. QSO 
clustering appears to be a realistic way to map large scale 
structure at high redshift, although the scales we are deal- 
ing with here are also difficult to study with existing samples 
(see, e.g.. Groom & Shanks 1996). However, analyses of the 
spatial structure of deep radio surveys (at a redshift z ~ 1) 
do appear promising (Gress et al 1996, Loan, Wall & Lahav, 
1997). 

Unless otherwise stated, the above cosmological param- 
eters will be used throughout and X-ray fluxes will be re- 
ferred to the 2-10 keV bandpass. 

1.2 Deep versus wide-area surveys 

The next generation of X-ray instruments to be operating 
in the next decade are mostly based on X-ray imaging with 
grazing-incidence X-ray telescopes working up to energies 



~ lOkeV. This is indeed necessary to unveil the origin of 
the XRB, since it should not be forgotten that most of the 
XRB energy density resides at 30keV, whereas most of our 
knowledge on the source content is so far limited to energies 
below 3keV where only a few per cent of the total energy 
budget is contained. It is then essential to carry out imaging 
surveys at as high an energy as possible and to identify the 
sources dominating the source counts. 

X-ray imaging telescopes are well-suited to medium and 
deep surveys. Their limited effective collecting-area calls for 
long integration times to reduce photon counting noise, re- 
sulting in very deep images, eventually down to the con- 
fusion limit. We discuss some of these missions and their 
relation to the present work in the next subsection. The 
source counts (the so-called log N - log S relation) are then 
determined down to very faint fluxes in a small solid an- 
gle and therefore with limited precision. If exposures are 
long enough so that confusion noise dominates over photon 
counting noise, a fluctuation analysis provides an extension 
of the log N - log S curve for almost another decade in flux 
downwards. 

An alternative way of using an X-ray telescope with a 
suitable imaging detector is by doing a shallow all-sky sur- 
vey. This was done by ROSAT at soft X-ray energies (Voges 
1993) and wiU be done by ABRIXAS up to 10 keV (Friedrich 
et al 1996). Typical exposure times are then short and pho- 
ton counting noise is the limiting factor. Such an experi- 
ment is able to produce a map with accurate positions of the 
brightest sources over the whole sky. One can then extract 
much information about the cosmographical distribution of 
the sources and the local LSS of the Universe. 

The method we propose here requires the whole sky to 
be surveyed (to achieve the best statistics) to measure the 
power spectrum of the density fluctuations (PS) , with neg- 
ligible photon counting noise. In order to avoid the strongest 
structures associated with the galaxy, only photon energies 
above 2 keV will be considered (see discussion in section 
4). A diffuse galactic component, amounting to < 10 per 
cent, associated to the galaxy has been detected (Warwick, 
Pye & Fabian 1980, Iwan et al 1982), but it is expected 
to be smooth on scales of a degree and therefore will not 
contribute to the fluctutaions. We therefore focus on non- 
imaging instruments in order to obtain a large effective col- 
lecting area, such that in a typical exposure of about 100 s, 
confusion noise (i.e., the fluctuations in the sky brightness 
caused by the presence of the sources) exceeds photon count- 
ing noise. The fleld of view cannot be collimated to much less 
than 1 deg^ to maintain the required low photon counting 
noise level. The most clear example of such an experiment 
was the HEAO-1 A 2 experiment, and, as an instrument, the 
Ginga Large- Area Proportional Gounter (LAC), which, un- 
fortunately for the purposes of this paper, did not carry out 
an all-sky survey. 

For a homogeneous distribution of sources, the sky 
brightness fluctuations on scales of a few degrees will 
be dominated by relatively bright (~ 10~^^ ergcm"^ s~^) 
nearby sources (see section 2 for details) . However, clustering 
of sources, if strong enough, can be visible out to more dis- 
tant redshifts. The reason is that if several or many distant 
faint X-ray sources cluster within the scale of a field-of-view, 
they will produce a large enough signal in the distribution of 
X-ray sky intensities. As it will be shown later, the imprint of 
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inhomogeneities in the distribution of X-ray sources on the 
excess fluctuations is weighted by the redshift-dependence of 
their X-ray volume emissivity. This means that fluctuations 
in excess of confusion noise produced by the bright, just 
unresolved, sources will be most detectable at the redshifts 
where the bulk of the XRB originates. 

It should also be emphasized, and this has been proven 
by similar analyses on existing data, that besides having a 
negligible photon counting noise, what determines the sensi- 
tivity in the measurement of excess fluctuations is the num- 
ber of independent observations - as it will be shown later 
i^—f-^^ oc N~^J'^ , Nobs being the number of sky positions 
where the XRB intensity has been measured. All-sky cover- 
age is then essential. 

1.3 X-ray missions relevant to the present worlc 

There are a number of previous, existing and planned mis- 
sions which have made or are expected to make decisive 
steps towards the goal that is pursued here. Among them, 
the HEAO-1 A2 all-sky experiment has proven to be the 
most useful for cosmological work (Boldt 1987), since it pro- 
vided all-sky coverage with small photon counting noise. X- 
ray HEAO-1 A2 maps have been used to measure the XRB 
dipole (Shafer & Fabian 1983, Shafer 1983, Lahav, Piran & 
Treyer 1997), maybe also higher order multipoles (Lahav, 
priv comm), the search for the Great Attractor (Jahoda & 
Mushotzky 1989), cosmography of voids (Mushotzky & Ja- 
hoda 1992) and superclusters (Persic et al 1990) and many 
other cosmologically relevant issues. Many things have been 
learned from that experiment, and in particular that an ab- 
solute determination of sky brightness requires a combina- 
tion of two different flelds-of-view. It will be shown here that 
when other data become available, the HEAO-1 A2 obser- 
vations will enable a decisive step forward to be made in the 
measurement of the PS at intermediate redshift. 

As mentioned before, the Gmga LAC, with a field-of 
view of ~ 1° X 2° and a larger effective area, could have pro- 
vided a very accurate measurement of the PS at high redshift 
if it had carried out an all-sky survey. Also, it was rather 
unfortunate for the present purposes that the Gmga LAC 
had all collimators with the same angular size and there- 
fore the non-cosmic XRB had to be modelled (as oposed 
to subtracted). Nevertheless, even with the data available, 
very interesting constraints on LSS were found (Carrera et 
al 1991, Carrera et al 1993). 

In this paper we propose to measure or constrain the 
excess fluctuations by improving the precision of the con- 
fusion noise produced by relatively nearby bright sources. 
That means that the source counts down to fluxes < 
10~^^ ergcm"^ s~^ need to be obtained with high accuracy. 
This relation is only known today from fluctuation anal- 
yses of HEAO-1 A2 (Shafer 1983) and Gmga LAC data 
(Hayashida 1989, Butcher et al 1997). More recently ASCA 
surveys (Inoue et al 1996, Georgantopoulos et al 1997) have 
measured the source counts at fluxes ~ 10"^'' ergcm"^ s~^. 
However, these results are based on small solid angles and 
therefore the normalisation of the source counts is uncertain. 

A forthcoming mission that will deflne the log N - log 
S relation at lO"^'^ ergcm"'^ s~^ fluxes with the highest ac- 
curacy is ABRIXAS. This happens because ABRIXAS wiU 
detect all sources in the sky brighter than this flux. A dis- 



cussion of the accuracy in the modelling of confusion noise 
enabled by ABRIXAS is presented in the next Section. 

The next most important mission for the proposed ex- 
periment is XMM. XMM will not only deflne the log N - 
log S at fainter fluxes, but, most importantly, will flnd the 
X-ray spectrum of the sources that contribute both to the 
confusion noise and the excess fluctuations down to very 
faint fluxes. XMM will also be measuring or constraining the 
power spectrum of the fluctuations at comoving wavevectors 
kc ^ 0.1 — 1 /iMpc~^, thus complementing the observations 
at larger scales. 

Finally, AXAF may resolve the whole XRB at ~ 1 keV. 
The superb X-ray angular resolution combined with large 
optical telescopes will provide a direct insight into the evo- 
lution of the X-ray volume emissivity of the sources produc- 
ing virtually all the XRB, which is one of the key inputs to 
model excess fluctuations in terms of source clustering. 



1.4 Organisation of the paper 

In the next section we parametrise the noise components 
and sensitivities relevant to the measurement of excess fluc- 
tuations by a coUimated field-of-view proportional counter. 
Photon counting noise is estimated in terms of cosmic and 
non-cosmic backgrounds, based on Ginga LAC and RXTE 
observations. Confusion noise is modelled according to our 
best (rather inaccurate) knowledge of the source counts at 
10"^'^ ergcm"^ s~^. Particular attention is paid to the is- 
sue of the uncertainties and biases that source variability 
might introduce in that modelling. We show that {^^^ ~ 
2 X 10^^ might be reachable with the HEAO-1 A2 maps 
when ABRIXAS source counts become available. 

Section 3 presents the relation between excess fluctua- 
tions and PS. Assuming that the 2-10 keV redshift evolution 
of the volume emissivity is similar to the one associated to 
the resolved component of the soft XRB (mostly contributed 
by QSOs and NLXGs), we demonstrate that sensitive mea- 
surements of the PS at 2 ~ 1 — 2 on comoving wavevectors 
kc ~ 0.01 — 0.1 /iMpc~^ can be achieved with a few degree 
collimator, such as HEAO-1 A2. This corresponds to a PS 
signal below that expected from Cold Dark Matter mod- 
els, which could therefore be accurately tested when reliable 
information on the 2-10 keV X-ray volume emissivity (pre- 
sumably from AXAF optical identiflcation of deep flelds) 
becomes available. 

Section 4 presents a more detailed study on what could 
be achieved with a new X-ray mission surveying the whole 
sky with a ~ 0.5 — 2deg^ beam. Using expected future 
information on the X-ray volume emissivity, such an ex- 
periment would be able to measure the PS at intermedi- 
ate redshift with great accuracy at comoving wavevectors 
kc ^ 0.01 — 0.1/iMpc~^. Specifically, constraints on the evo- 
lutionary parameters of the PS (and in particular go) at the 
involved redshifts could be found. Such a mission would pro- 
vide other scientific results which are also outlined in that 
section. 

In Section 5 we summarize our results. 
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2 SENSITIVITY OF XRB OBSERVATIONS TO 
EXCESS FLUCTUATIONS 

For a uniform distribution of X-ray sources, the distribution 
of spatial fluctuations in the XRB, when observed through 
a given beam, can be predicted if the source counts and the 
noise components (particularly photon counting noise) are 
known. This is well documented (Scheuer, 1957, 1974; Con- 
don 1974, Shafer 1983) and has been applied many times 
with success to derive X-ray source counts from XRB fluc- 
tuations (Shafer 1983, Hamilton & Helfand 1987, Barcons 
& Fabian 1990, Hasinger et al 1993, Barcons et al 1994, 
Butcher et al 1997). Under the assumption of a uniform dis- 
tribution of sources in the sky and of confusion noise (i.e., 
the noise coming from the presence or absence of sources) 
dominating over photon counting noise, the distribution of 
XRB fluctuations (the so-called P(D) curve) is only sensi- 
tive down to a flux where there is about one source per 
beam, and its width (standard deviation) is dominated by 
the sources of which there are a few (this number depends on 
the slope of the source counts, see below) per beam. When 
photon noise is important, the total intrinsic dispersion of 
the intensity distribution is then 



(-) 

\ I J intrinsic 



Ixi 



(1) 



where CTc is the confusion noise, Oph is the photon counting 
noise and Ixrb is the XRB intensity. As a practical issue, 
we should emphasize that an overall fit to the P(D) curve 
is mostly sensitive to its second moment, but it has the ad- 
vantage over the variance that it is not dominated by the 
brightest unresolved sources (i.e., the tail of the P(D)), but 
by the sources where there are a few per beam. 

Now, when the sources are not uniformly distributed, 
but are clustered in the sky, the situation is different. The 
predicted P(D) distribution needs then all of the n-point 
correlation functions, for n at least as large as the total av- 
erage number of sources per beam (Barcons 1992). There are 
indeed models for the n-point correlation function of objects 
(e.g., a simple gaussian model where all correlation functions 
beyond n = 2 are zero, or a model where there is a ran- 
dom distribution of clusters of sources, all of them with the 
same average profile which determines the n-point correla- 
tion functions) , but observationally there is little knowledge 
of them beyond n = 3. 

The contribution of clustered sources to the shape of 
the P(D) does not affect only the sources brighter than the 
one-source-per-beam level, but all sources equally. Clusters 
of very faint sources might produce significant dispersion 
in the XRB intensity. As will be shown later, the effect of 
clustering is weighted by the X-ray volume emissivity as a 
function of redshift. 

Fortunately, in most situations the P(D) shape is domi- 
nated by the confusion noise of relatively bright sources, and 
the effect of clustering is only a small correction. In those 
cases it has been customary to parametrise the effect of clus- 
tering in terms of 'excess fluctuations', i.e., a small quantity 
to be quadratically added to the intrinsic dispersion 
of the P(D) (eqn. 1). In what follows we shall assume this 
approach. 

The key point of this paper is based on the fact that the 
measurement of the variance of an approximately gaussian 



distribution (we can use that approximation for P(D) for 
this particular purpose) is distributed as . Therefore the 
2(7 uncertainty with which the dispersion of the P(D) can 
be measured is {2/ NobsY^'^ times the dispersion itself, where 
Nobs is the number of independent observations of the XRB 
intensity. That means that the 2(j sensitivity at which excess 
fluctuations could be measured is 



2 (ag-Hcr;^)^ 

Nobs IxRB 



(2) 



If we have a beam with solid angle Hdeg^, and the whole 
high galactic latitude sky (| b |> 20°) is used, then this 
number is ~ O.OIQ^^^ (^) , which is why all-sky 

\ ^ / intrinsic 

coverage is essential. 

It should be emphasized that eqn. (2) shows the max- 
imum precision that can be reached and it requires knowl- 
edge of the intrinsic dispersion (eqn. 1) to better than this 
flgure. In what follows we estimate these values under gen- 
eral grounds. 

We assume a proportional counter with coUimated field 
of view of solid angle Odeg^, effective area 10*^4cm^ with 
an energy bandpass from ei to 62- For most purposes we will 
use £1=2 keV and £2 = 10 keV. This instrument is assumed 
to scan the whole sky in 6 months, so the typical exposure 
time will be of the order oi t = lOOtioo^^^''^ s. 



2.1 Photon Counting Noise 

Proportional counters detect events from the cosmic XRB 
as well as from particles crossing the detector. These have 
different energy spectra: the XRB has an energy spectrum 
Ixrb oc e""'*, which should be folded through the appro- 
priate instrumental response (assumed as a constant effec- 
tive area here) to be converted to counts, while the particle 
background is often well approximated by a constant energy 
dependence in the number of counts detected. If we write the 
total count-rate as the sum of these two terms 



Cxot = CxRB + CdET 

then 

CxRB = ain44*(-0.4,ei,e2) 
and 

Cdet = a2/l4*(l, ei, £2) 



(3) 
(4) 
(5) 



where we take into account bandpass corrections through 
the function 



*(Aei,e2) = 



lof - 2" 



(6) 



The constants ai and 02 are more or less universal for the 
same type of orbit. We have estimated their values based 
on the Gmga LAC and the RXTE PCA. For Gmga LAC 
{A4 = 0.4, t7 = 2, £1 = 4, £2 = 12), CxRB ~ 8cts-^ and 
Cxot ~ 14cts~^ (see, e.g., Carrera et al 1993), while for 
RXTE {Ai = 0.14, J7 = 1, £1 = 2, £2 = 10), Cxrb ~ 
2.5 cts^^ and Grot ~ 5cts~^ (K. Jahoda, private commu- 
nication). From both instruments we find consistent values 
around ai ~ 02 « 15 — 17cts~'^. 

The photon counting noise contribution to the disper- 
sion of the P(D) curve can then be estimated as 
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O'ph 
V IxRl 



0.024 {Aitmoy^ x 



+ 



3'(l,ei,e2) 



(7) 



The accuracy of this value depends crucially on how 
stable the particle background is around the orbit and on 
whether it can be subtracted rather than modelled (this re- 
quires collimators of various sizes). An equatorial orbit min- 
imizes variations due to the particle background. 



2.2 Confusion Noise 

In what follows we assume that the source counts dominat- 
ing the confusion noise follow a euclidean power law 



— = ;fn(7-i)— f— 

dS So V So 



(8) 



where 7 = 2.5, 5*0 is a reference flux arbitrarily chosen 
to be 10"^* erg cm~^ s~^ (2-10 keV) and the normalisation 
K = 300 7^300 sources deg~^ is the number of sources per 
square degree brighter than So- Both the slope and the nor- 
malisation of the source counts are consistent with the fluc- 
tuation analysis carried out with the Ginga LAC (Butcher 
et al 1997) and also match the deeper ASCA surveys (In- 
oue et al 1996, Georgantopoulos et al 1997), all of them 
referred to the 2-10 keV band. The spectrum of the Ginga 
LAC fluctuations also showed that the sources responsible 
for the confusion noise on 1 deg^ scales have a power-law 
spectrum with energy spectral index 0.7 and negligible ab- 
sorption (Butcher et al 1997). 

Confusion noise can be estimated following Condon 
(1974). If beams with an intensity more than V<j above the 
mean are removed since they will be identified as sources, 
the variance of the remaining map can be estimated itera- 
tively and the confusion noise (defined as the fiux equivalent 
to a Icr signal in the intensity histogram) is 



So KQ. 



7- 



1, 



,3-7 



(9) 



which for a typical value F « 5 and euclidean counts is 

(Tc ~ 1.6 X 10"^^ergcm"^s"^ (ifsoo^^)^ • (10) 

This confusion noise produces a contribution to the intrinsic 
dispersion of the P(D) curve 



IxRB 



Q.nKiooO. 



»(-0.7,ei,£2) 
*(-0.4,ei,£2) 



(11) 



which will have to be added in quadrature to eqn. (7) to find 
the total intrinsic dispersion of the P(D) curve. 

A very important point is to what degree of accuracy 
this term can be estimated, since it is likely to be the main 
limiting factor to measurement of excess fluctuations. It is 
indeed important to have an all sky survey (as ABRIXAS 
will provide) in order to minimize the statistical errors on 
this quantity. Over 8500 7^300 sources are expected at high 
galactic latitudes down to a flux of 10~^^ erg cm~^ s~^ which 
would mean an accuracy of the order of 1 per cent in K and 
slightly better for CTc- 

However, there is the issue of the source variability. 
AGN, which are supposed to be the dominant class of source 
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Figure 1. Relative shift in the estimates of confusion noise from 
an all-sky survey as a function of the slope of the differential 
source counts. Filled dots correspond to a factor of 2 variability 
and triangles to a factor of 5 variability (these have been slightly 
shifted to higher values of 7 to avoid confusion between error 
bars). The error bars represent the standard dispersion due to 
statistics and variability. 



at a flux 10~^^ ergcm~2 are known to vary substantially 
on all timescales. Assuming that sources vary independent 
of each other and with a similar (flux-independent) ampli- 
tude, the measured source counts in a flux limited sample 
are slightly different from the average source counts (which 
are relevant to the P(D)). The main reason is that the steep- 
ness of the source counts makes a fraction of the sources of 
average flux below the detection threshold contribute to the 
source counts above the detection threshold. Indeed some 
of the sources of average flux above the detection threshold 
would also be undetectable, but these are less numerous. The 
net result is that the source counts in a flux limited sample 
are an overestimate of the average source counts, and it is 
this last one which matters for the confusion noise. 

In order to quantify the variability effect on the confu- 
sion noise we have simulated fluxes of sources, whose aver- 
age values are drawn from the source counts given by eqn. 
8 down to 3 X 10~^'^ergcm~^ s~^. The number of simulated 
sources is set to cover the whole high galactic latitude sky 
in order to mimic the ABRIXAS survey. The fluxes are then 
allowed to vary randomly within a factor of several. That 
produces a different list of source fluxes which is then cut 
at 10~^^ ergcm"^ s~^. The source counts are then fltted via 
maximum likelihood to a single power law (as in eqn. 8), 
resulting in an estimate of the confusion noise given in eqn. 
9. 

Fig. 1 displays the relative variation between fitted and 
expected confusion noise as a function of input slope 7 for 
^^^300 = 1 when all sources vary within a factor of 2 and a 
factor of 5. The amount of overestimation of the confusion 
noise is ~ 1 per cent for euclidean counts if sources vary 
within a factor of 2 and, as expected, it grows with 7. Since 
not all of the sources vary (there will be some contribution in 
the source counts from clusters of galaxies, galaxies, etc.) the 
error is likely to be smaller. For euclidean counts the error is 
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variances as small as 2 x 10 ^ could be measured even with 
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Figure 2. Minimum detectable excess variance for all-sky obser- 
vations and different beam sizes 



therefore small enough to reach the maximum accuracy in 
the excess fluctuations, assuming source variability within 
a factor of 2. The situation is clearly worse if sources vary 
within a factor of 5, where for euclidean counts the variation 
in the confusion noise will be of the order of 5 per cent. There 
are indications (R.S. Warwick, private communication) that 
source variability could be large on average (close to a factor 
of 5) from the comparison of Rosat All-Sky Survey data 
with later pointed observations. However, since variability 
produces a systematic effect on the confusion noise, different 
observations carried out by ABRIXAS (e.g., observations 
taken 6 months apart) of the same sources could be used to 
quantify this effect and to correct the source counts. 



2.3 Sensitivity to excess fluctuations 

Eqn. 1, when combined with eqn. 11, shows that the min- 
imum intrinsic dispersion in the P(D) curve for a ~ 1 deg 
set of observations is around 10 per cent. The accuracy with 
which any excess variance could be measured on top of this 
intrinsic value depends on both the precision with which this 
dominating term can be modelled. If the intrinsic dispersion 
can be modelled to better than 1 per cent (which is close to 
the maximum statistical accuracy in the absence of source 
variability), then excess variances as low as ~ 10"'^ could be 
detected. 

More precisely, and taking into account the estimated 
values for the confusion and photon counting noise from the 
previous subsections and assuming that source variability is 
not going to dominate the precision with which confusion 
noise can be estimated, the excess fluctuations that could 
be eventually measured by an all-sky observation with a 
collimated field-of-view proportional counter are shown in 
Fig. 2. This assumes a 2-10 keV bandpass. 

Clearly, when the product ^4tioo is larger than a few 
(the precise value depending on the beam f2), the intrin- 
sic dispersion of P(D) is dominated by confusion noise as 
opposed to the photon counting noise which dominates at 
smaller values. When P(D) is confusion dominated, excess 



relatively large beams. For Q, 
would be within reach. 



1, values close to 1 x 10 



3 MEASURING THE POWER SPECTRUM OF 
THE FLUCTUATIONS AT HIGH REDSHIFT 

The excess variance can be related to the power spectrum of 
the density fluctuations in the following way (see Barcons & 
Fabian 1988, Carrera, Fabian & Barcons 1997 for details): 



(12) 



d'qG^q) {^^y V{z-q/dA{z)) 
where the XRB intensity Ixrb is 



IxRB ~ 



•ff 



47r Ho 



dz{l + z)-''{l + 2qoz)-^j{z). (13) 



fleff = 3.046 X 10~*f2, j{z) is the X-ray volume emissivity 
at redshift z (with the appropriate K-correction), dA is the 
angular distance 



dA(z) 



c [zqo + (qo-l)i-l + Vl + 2qoz)] 



(14) 



G{q) is the 2D Fourier transform of the beam function and 
'P{z; k) is the power spectrum of the fluctuations (fc = (1 -|- 
z)kc is the physical wavevector) , which is related to the 
source 2-point correlation function £^{z; r) by 







-Yv{z-k) 



(27r) 



d^re ^'^'^^{z;r), 



(15) 



and finally 



Fiz- go) = (1 + ^)-'(l + 2qoz)-^f(z)/dliz). (16) 

The above equations provide the basic link between excess 
fluctuations, evolution in the X-ray volume emissivity (the 
actual normalization is irrelevant since it cancels out) and 
the power spectrum. Although the PS enters in the expres- 
sion of the excess fluctuations in a convoluted way, there are 
a couple of simplifications that provide an almost one-to-one 
relation between power spectrum and excess fluctuations. 

In what follows a circular beam with gaussian proflle 
(dispersion angle ^o) will be assumed. 

The first simplification to realise is that only wavevec- 
_ i_ 

tors k ~ Q.^j/'j: dA{z)~^ axe relevant. The effective filtering 
function 2-KqCP'{q) is shown in Fig. 3 for a Ideg FWHM 
beam (i.e., 6o — 1/2.354 deg). It can be seen that only val- 
ues of the 2D wavevector q around the maximum qmax = 
2~2g~^ 100r2~2 deg~^ will contribute to the excess fluc- 
tuations. That implies that at every redshift z, only the 3D 
wavevectors k ~ qmax / dA{z) will contribute substantially to 
the excess fluctuations. For beamsizes of the order of a de- 
gree, and signiflcant redshifts 2: ~ 1, the 3D wavevectors to 
which it is maximally sensitive are k ~ 0.01 — 0.1 /iMpc~^, 
which is near the expected maximum of the PS. Then, it 
is safe to assume that the PS does not vary much over the 
range of relevant wavevectors, which results in the following 
expression for the excess fluctuations 
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Figure 3. Effective beam filtering function for a 1° FWHM circu- 
lar beam with gaussian profile as a function of the 2D wavevector. 
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Figure 4. The weight W{k, z) defined in eqn. 19 as a function of 
redshift for representative values of k 
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where 



W{z) 



2Ho (1 + zy'il + 2qoz)-2f(z)/d\{z) 
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(17) 



(18) 



The next simplification is hinted at by our (limited) 
knowledge of the redshift evolution of the X-ray volume 
emissivity j{z). It is hoped that after AXAF and XMM are 
launched, and deep surveys have been carried out, optical 
identification work (especially for the AXAF sources whose 
positions will be determined with very good accuracy) will 
be able to reveal the volume emissivity j{z) as a function of 
redshift. So far, at soft energies it appears that AGNs have 
a steeply rising volume emissivity {j{z) oc K{z){l + z)^'^'', 
K{z) being the K-correction and p ~ 3) up to Zc — 1.7, but 
beyond that redshift everything is consistent with no evolu- 
tion (Boyle et al 1994, Page et al 1996). The NLXGs that 
appear to be more numerous at faint fluxes, might also show 
a similar behaviour (Griffiths et al 1996). To our present 
knowledge then, j{z) is a steeply rising function up to some 
redshift Zc and then it flattens out to a constant comoving 
volume emissivity. 

In practice that means that the influence of the power 
spectrum (i.e. inhomogeneities in the distribution of sources) 
on the excess fluctuations is heavily weighted towards Zc for 
wavevectors k < 0.1/iMpc~ . In order to illustrate this fact, 
we define the function W{k, z) 



^ dk dzW{k,z)P{z;k) 



(19) 



Fig. 4 shows W{k, z) for different values of k assuming p = 3, 
an energy spectral index a = 1 for the sources to com- 
pute the K-correction (the actual energy spectral index is 
likely to be smaller and therefore the function W{k,z) will 
be much more peaked towards Zc) and Zc = 1.7. What is 
seen there is that the wavevectors that dominate are around 
k ~ 0.1/iMpc~^ and that within these, it is the redshift 



at which the volume emissivity peaks that is most heavily 
weighted. There is of course some contribution from smaller 
wavevectors at lower redshifts, but since the peak in the 
power spectrum is expected to be found around these shorter 
wavevectors, not much contamination from large-scale local 
structure is to be expected. In fact, this large-scale power 
could be removed by 'flat-fielding' the all-sky maps with 
a large-scale smoothed version of the same maps, leaving 
scales of ~ 0.1/iMpc~^ unaffected. This would have the 
additional advantage of removing any residual galactic large- 
scale structure. It is then concluded that we could be mea- 
suring the PS at a redshift beyond the deepest available 
galaxy surveys. 

To estimate the sensitivity in the measurement of the 
PS with X-ray observations, we assume that all of the XRB 
comes from a redshift bin Az = 2 around Zc = 1.7 and 
approximate the integrals in redshift as the central value of 
the integrand times Az. We then find 



f — V - — 



\ I J 



AV 



T{zc; ku), 



(20) 



where ko = 2 ^9^ ^dA{zc) ^ and AV^ is the volume sampled 
by a beam 



AV = n,ffdA{zcf cHo\l + z)-^{l + 2qoz)-iA2 



(21) 



Fig. 5 shows the maximum sensitivity, in terms of the 
power spectrum as a function of comoving wavevector, for 
different beam sizes, according to the estimates from Sec- 
tion 2. A factor of 2 reduction in the number of indepen- 
dent measurements has been included in order to account for 
the fact that neighbouring observations will not be indepen- 
dent. For the HEAO-1 A2 points (triangles) AaIioq = 0.1. 
The dashed lines represent the expected maximum sensitivi- 
ties for all-sky observations of the XRB with different beam 
sizes (ranging from lOdeg^ to 0.1 deg^ and various values 
of AiXioo- The filled dots represent beams of 2, 1 and 0.5 
deg^ for AitwQ = 1, as in the instrument whose concept is 
presented in the next section. From this figure it is clearly 
seen that the best sensitivity near the expected peak of the 
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Figure 5. Power spectrum of density fluctuations as a function of comoving wavevector all computed at a redshift z = 1.7. The solid line 
represents the prediction from local observations evolved with go = 0-5- The dashed curves are the 2(r sensitivities from all-sky surveys 
carried out with various beam sizes and for the quoted values of the parameter AntiQQ. The triangles show the maximum sensitivity that 
could be reached with HEAO-1 A2 data (for the 3° X 3° and 1.5° X 3° collimators from left to right) and the filled points show the 
sensitivity of the instrument whose concept is presented in Section 4.1, with beam sizes of 2, 1 and 0.5 dcg (left to right). The dotted 
line shows the sensitivity of XMM observations for 2 years as explained in the text. 



PS is achieved by 0.5-1 deg^ beams and that to avoid the 
results being dominated by photon counting noise (i.e., the 
flattening of the dashed curves towards high wavevectors) 
at those angles, a value of Aitioo close to 1 is needed. 

In the same figure we also show the expected sensitivity 
reached at larger wavevectors with XMM observations over 
2 years. That has been computed assuming a log N -log 
S as observed at soft X-ray energies with 7 = 2 but with 
a normalisation two times larger, as it seems to apply to 
the 2-10 keV passband (Georgantopoulos et al 1997). We 
have assumed about 500 useful observations with an average 
exposure time of 20 ks during that period. 

We also show for comparison a Cold Dark Matter power 
spectrum (see Peacock 1997 and in particular Fig. 6 of that 
paper) linearly evolved to redshift 1.7 (with go = 0.5). A 
shape parameter F* = 0.25 has been assumed, claimed by 
Peacock & Dodds (1994) to fit the shape of the local LSS 
very well and the normalisation has been chosen accordingly. 
The conclusion is that when the confusion noise from sources 
brighter than ~ 10~^^ ergcm"^ s~^ can be accurately mod- 
elled, the all-sky HEAO-1 A2 observations might be sensi- 
tive enough to measure the PS at intermediate redshift (that 
is in the absence of other systematics). With a smaller beam 
(close to 0.5 — Ideg^) and larger area the accuracy in the 
measurement could be of a few per cent. 



4 MEASURING THE XRB FLUCTUATIONS 
ON 1 DEGREE SCALE 

4.1 A mission dedicated to the XRB 

Throughout this paper we have illustrated what results 
would be obtained by using a beamsize close to 1°. In what 
follows we want to show how the simplest possible instru- 
ment (i.e. a collimated field-of-view proportional counter) 
could provide a very significant cosmological result. 

Indeed other instruments might provide in principle 
similarly valuable information with comparable sensitiv- 
ity in terms of excess variance. For values of the product 
^4^100 > 0.1, the measurement of the fluctuations could be 
sensitive enough to provide cosmologically relevant results. 
We believe that systematics are going do dominate the ul- 
timate sensitivity of such measurements. If unpredictable 
and slow gain variations in the instrument are the domi- 
nant source of the systematics, then the larger the effective 
area the better, especially when the measurement can be re- 
peated a few times. A detailed study would be required for 
each experiment where systematics need to be understood 
and kept to a minimum. We devote here special attention 
to proportional counters which are well understood and and 
have proven stable over extended periods. 

Although the smallest detectable excess fluctuations 
(■^V for a 1 deg^ expriments would be within a factor 
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Figure 6. Contour levels of constant power spectrum (in 
Mpc'^ units) achievable at redshift Zc = 1.7. The solid lines 
represent the 2-10 keV bandpass and the dotted lines the 1-5 keV 
bandpass. The 2 HEAO-1 A2 collimator sizes (4.5 and 9 deg^) 
have been labelled "A2" , while the instrument proposed here has 
been labelled "X" (collimator sizes 1 and 2 deg^). The ABRIXAS 
point, assuming an area of 25cm'^ and a total integration time of 
4000 s is shown as 'AB'. 

of 2 of what can be done with HEAO-1 A2 (see Fig. 2), the 
sensitivity in the measurement of the power spectrum would 
be increased by a much larger factor, since a much smaller 
volume AT/ would be sampled by a single beam (see eqn. 
20). To further illustrate this point we show in Fig. 6 what 
such an instrument would be able to do in terms of mea- 
suring the power spectrum, as a function of the parameter 
Aitioo and beam Q.. Fig. 6 also shows the expected values 
for a 1-5 keV bandpass which would collect more counts, but 
due to the fact that the sources that dominate the confu- 
sion noise have a steeper spectrum than the XRB, the P(D) 
curves would be noisier and the experiment less sensitive. 

We propose an instrument of effective area Ai — 1, with 
two collimator sizes (1 and 2 deg^) so the contribution from 
cosmic and detector backgrounds could be well separated. It 
might be actually very interesting that the collimators are 
elongated (e.g., 0.5° x 2° and 1° x 2°) in which case there 
will be some information on the power spectrum down to 
smaller scales. With such a large effective area, the X-ray 
brightness of the sky at any point could be determined with 
an accuracy better than 2 per cent. The power spectrum at 
kc ~ 0.01 — 0.1 /iMpc~^ could then be measured, in princi- 
ple, down to 0.1/i~''Mpc^ or better {2a) which would guar- 
antee not only a detection (if the PS at intermediate redshift 
is not largely overestimated by the above parametrisation) 
but also an accurate measurement. In what follows we try 
to highlight the need for such an instrument and also list a 
few complementary studies that would benefit considerably 
from such observations. 

First of all. Figs. 5 and 6 show the minimum detectable 
power spectra, assuming that the intrinsic width of the P(D) 
can be determined to its absolute statistical precission. In- 
deed, it might be that source variability is found to have a 
larger effect or that photon counting noise cannot be mod- 



(h Mpc 

Figure 7. Dimensionless power spectrum estimated at different 
values of qg. The solid curve corresponds to a linearly evolved 
go = 0.5 local PS (see text). The dashed line shows the same 
thing for go = 0.1. The filled and empty dots show the expected 
2(T sensitivities from the observations suggested in Section 4.1 and 
for go = 0.5 and go = 0.1 respectively. 

elled to 1 per cent accuracy. If one or both of these effects 
cause a 5 per cent uncertainty in the intrinsic P(D) noise, 
then with the HEAO-1 A2 experiment (degrading its sensi- 
tivity in the power spectrum by almost an order of magni- 
tude) the cosmic signal in the PS would be missed. A 5 per 
cent accuracy in the intrinsic P(D) width with the proposed 
experiment would nevertheless allow a power spectrum as 
small as ~ 2h~^ Mpc^ (still below the prediction) to be de- 
tected at 2a. Indeed, the fact that such an instrument would 
be working close to the peak of the power spectrum helps. 

To further illustrate the capabilities of this observation, 
we explored the possibility that a sensitive measurement of 
the PS at intermediate redshift could constrain the evolu- 
tionary parameters of the PS. We have taken the z = Cold 
Dark Matter spectrum with constant F* (see Section 3 and 
Peacock 1997) and evolved it linearly (which is appropriate 
to the scales under consideration) to 2; = 1.7 with different 
values of go. The dimensionless power spectrum {4TvklV{kc), 
Fig. 7) does not depend on the Hubble constant h, but does 
depend on go. The expected 2a accuracy for the dimension- 
less PS is between one and two orders of magnitude smaller 
than the predictions. The difference between the expected 
dimensionless power spectrum for qo — 0.1 and qo = 0.5 is 
only a factor of 4 at that redshift and therefore go could be 
measured in principle. 

Such a mission makes no severe technical demands re- 
quiring only a continuous scan of the whole sky over 2 years 
with the most stable proportional counters and 1-axis sta- 
bilisation. In order to keep the contribution to the intrinsic 
dispersion of the P(D) curve well determined, it is necessary 
to have the detector background as stable as possible. In 
this respect, the mission concept we propose here would be 
better suited if it is launched into an equatorial orbit (as for 
BeppoSAX). 
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4.2 Other ways to measure Large-Scale Structure 

There are many other scientific goals that could be achieved 
while such an instrument is performing its main task. Within 
the subject of the LSS, the multipoles of the XRB should 
be mentioned first. In fact at least the dipole needs to be 
well measured and compared with the dipole in the distri- 
bution of the X-ray sources to estimate the bias parameter. 
The earliest attempts to measure LSS with the use of the 
XRB looked also for 12-hr and 24-hr effects which could be 
attributed to the presence of a large lump of X-ray emit- 
ting matter at different distances (Warwick, Pye & Fabian 
1980). HEAO-1 A2 provided the possibility of a measure- 
ment of the dipole of the XRB, whose direction is in gen- 
eral agreement (within large errors) with the CMB dipole 
(Shafer 1983, Shafer & Fabian 1983). The amplitude of the 
XRB dipole is expected to be larger than the CMB one pri- 
marily because of aberration ((■%^)jj. = (3 + 0)7, Q being 
the energy spectral index which is a ~ 0.4 for the XRB and 
a = — 2 for the CMB), but also because the mass overden- 
sity which is gravitationally pulling the Local Group is also 
expected to emit X-rays above the average (Warwick, Pye 
& Fabian 1980, Miyaji 1995, Miyaji & Boldt 1990 and Boldt 
1992). 

More recently Lahav, Piran & Treyer (1997) have pro- 
posed a formalism to measure multipoles in the XRB and 
to relate them to the power spectrum of the fluctuations, 
as is done with galaxy surveys. Of course the main prob- 
lem in measuring multipoles is not a sensitivity one, but 
confusion noise. The amplitude of the dipole is < 1 per 
cent and must be measured in maps where confusion noise 
is larger. Since all the multipoles are in fact variances of 
the sky maps weighted with appropriate spherical harmon- 
ics, the contribution of confusion noise is dominated by the 
brightest sources that have not been removed. Than makes 
the measurement of multipoles particularly difficult (Lahav, 
private communication). In addition there is the unknown 
contribution of the galaxy, even at high galactic latitudes 
which could contribute to the low-order (large-scale) multi- 
poles. 

There is however an advantage in the experiment pro- 
posed here with respect to the only previous one that car- 
ried out an all-sky survey {HEAO-1 A2), when combined 
with ABRIXAS. The 'shot noise' variance for all multipoles 
is (see Lahav, Piran & Treyer 1997) 

(|a,,„|^).„ = 5o^i^nl-l(^)'"' (22) 

where Smax is the flux above which all sources have 
been removed. If ABRIXAS can flnd the positions of all 
sources brighter than lO^^^ergcm^^ s~^, of surface density 
is 0.3 deg~^, then one out of every 3 beams will have to 
be excluded from the multipole analysis. This still provides 
enough data, with virtually no impact on the amplitude of 
the multipole signal, but with the 'shot noise' reduced by 
a factor of 5, when compared with the HEAO-1 A2 maps 
where the Piccinotti et al (1982) sources have been removed. 

The proposed mission will also provide a good measure- 
ment of the Autocorrelation Function. Since the instrument 
would scan the sky along great circles, measurements of the 
sky brightness at separations less than the collimator fleld- 
of-view will be taken. These can be used to measure the 



autocorrelation function of the XRB on scales smaller than 
the beam size. However, since the beams will strongly over- 
lap in adjacent measurements, the signal in the Autocorre- 
lation Function will be dominated by the brightest sources 
present (see Carrera et al 1993 and Martm-Mirones et al 
1991). At the very least, this will be useful to conflrm the 
source counts at high fluxes (as was done, for example, in 
the Ginga High Galactic Latitude Survey by Kondo 1991). 
Beyond the angular scale where the beams overlap, the au- 
tocorrelation function is expected to measure true source 
clustering and extension of the cosmic sources. If clusters 
of galaxies can be reliably removed from the all-sky maps 
(again using ABRIXAS), source clustering on scales larger 
than the beam size could also be measured. 

Other possibilities for such an instrument include the 
search for positive or negative signals around known struc- 
tures (Great Attractor, Superclusters, Voids, etc.), studies 
of the cross- correlation function between XRB intensity and 
galaxy and cluster catalogues (Jahoda et al 1991, 1992; La- 
hav et al 1993a, Miyaji et al 1994, Barcons et al 1995, Car- 
rera et al 1995, Soltan et al 1996) as well as with X-ray maps 
at softer energies, cross-correlations of the XRB with Cos- 
mic Microwave Background maps (Boughn & Jahoda 1993; 
Boughn, Crittenden & Turok, 1997, Kneissl et al 1997) pos- 
sible detection of excess skewness in the fluctuations (similar 
to what is done in counts-in-cells, Lahav et al 1993b) and 
many types of studies of the distribution of X-ray sources 
and diffuse emission from the Galaxy. 



5 CONCLUSIONS 

We have shown that all-sky degree-scale observations of 
the XRB carried out with enough sensitivity (effective area 
times exposure time in excess of 10^ cm^s) could reveal ex- 
cess fluctuations due to the clustering of distant sources 
down to levels (■%^) ~ 1 — 2 x 10~^. These excess fluctua- 
tions are likely to arise from high redshift sources and then 
the power spectrum of the fluctuations at those early epochs 
could be measured with sufficient sensitivity to detect the 
cosmic signal. 

The generic requirement for such a goal is that the in- 
trinsic width of the P(D) curve (caused by confusion noise 
and photon counting noise) is estimated to better than 1 per 
cent. Source variability (if within a factor of 2) already in- 
troduces a bias of 1 per cent for euclidean source counts if an 
all-sky survey of sources brighter than 10"^'^ ergcm"^ s~^ is 
used (as will be carried out by ABRIXAS). If variability is 
as large as a factor of 5, then the contribution of confusion 
noise to the intrinsic P(D) dispersion will suffer from a 5 per 
cent error. Besides this, photon counting noise needs to be 
modelled to better than 1 per cent to achieve the maximum 
sensitivity in the excess fluctuations. 

When the X-ray source counts down to that flux have 
been measured, the P(D) noise in the HEAO-1 A2 maps 
could be accurately modelled and the power spectrum of 
the density fluctuations at intermediate redshift might be 
measurable if the modelling of P(D) can be done accurately 
enough. 

A new X-ray mission consisting of a large effective area 
proportional counter with 1 and 2 deg^ collimators, which 
would scan the whole sky several times would provide a 
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much more secure approach. Such an experiment would have 
a sensitivity 10 times better in terms of the power spectrum 
than the HEAO-1 A2 experiment and would be able to de- 
tect a signal in the power spectrum over 10-100 times smaller 
than the predictions of the popular models. We believe that 
such an experiment constitutes the best chance to measure 
the power spectrum of the density fluctuations in the Uni- 
verse at a redshift z ~ 1 — 2 where the different cosmological 
scenarios give the most distinct predictions. 
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